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We compare disentanglement and decoherence rates within two-spin and three-spin entangled 
■ systems subjected to all possible combinations of local and collective pure dephasing noise combi- 

nations. In all cases, the bipartite entanglement decay rate is found to be greater than or equal to 
, the dephasing-decoherence rates and often significantly greater. This sharpens previous results for 

, two-spin systems [Phys. Rev. B 56, 165322 (2003)] and extends them to the three-spin context. 
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I. INTRODUCTION 



> 

■ One goal of nanoscience is to understand the behavior of quantum entanglement in states of multiparticle systems, 
such as spin arrays. In such contexts (e.g. [l[), spins are often called qubits, as we will refer to them here. Entanglement 
can be viewed as a resource that allows uniquely quantum phenomena, such as quantum state teleportation [2| using 
quantum memories, to be exhibited [1, Hj]. Entanglement can be affected by a number of noise influences, including 
both vacuum noise [f| and classical noise [f|. One important effect is that of the loss of entanglement, which arises 
in both these cases. Recently, for example, it has been noted that in initially entangled two-qubit systems a sudden 
loss of entanglement can occur, an effect known as "entanglement sudden death" 0,0]. It is also interesting to note 
that in certain situations the irreversible process of spontaneous emission may lead to the revival of entanglement 
Q • An experimental method for the probing of disentanglement in cavity QED has also been introduced Q . Recent 
analyses have analytically examined the relationship between dephasing decoherence and the reduction of bipartite 
entanglement in two-qubit systems under such noise [t| [l(| El EH , as well as in pairs of three- level systems, i.e., 
two-qutrit systems [l3|. Here, we extend the analytic consideration of this relationship to the case of three-qubit 
^-J '■ systems. 

^ , Decoherence results from unwanted interactions of a quantum system with its environment. When such decoherence 
V "~j ■ occurs, it is exhibited by the decay of off-diagonal elements of the density matrix describing the system [l4[ . This may 
be due to either local or collective interactions, acting either at the scale of subsystems or of the composite system; 
that is, environmental noise may act on some or all system constituents similarly or differently and on some not at all 
0, E3] • Ex amp les of environments capable of inducing such interactions include electromagnetic fields and thermal 
baths [ll.l9l.ll0j] . Here, we consider quantum state decoherence in systems of two or three qubits due to pure-dephasing 
noise. In particular, we address the relationship between degree of entanglement (concurrence) and state coherence 
under environmental noise by determining characteristic decay timescales at the single-qubit, double-qubit, and triple- 
qubit scales. We find decoherence timescales to bound corresponding disentanglement timescales from above in all 
cases. 

We begin by reviewing and sharpening the analysis begun in an earlier study of two-qubit systems, and then extend 
that approach to the case of three-qubit systems, in which more complex combinations of dephasing noise are relevant. 
This paper is accordingly organized as follows. In Sec. II, we revisit the analysis of two-qubit systems begun by Yu 
and Eberly in [Toj] that explicitly treats disentanglement in the operator-sum formalism under local dephasing-noise 
environments and explicitly address disentanglement under collective dephasing noise, and sharpen the conclusion of 
the analysis by showing that decoherence times bound disentanglement times from above in all cases. In Sec. Ill, we 
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extend the analysis of the effects of dephasing on coherence and bipartite entanglement to three-qubit systems in the 
more complex dephasing noise environments available there. We conclude with a summary of results in Section IV. 



II. TWO-QUBIT SYSTEM 



Before considering the effects of dephasing on three-qubit systems, we complete a previous analysis of two-qubit 
systems by Yu and Eberly [10( • There, a number of dephasing scenarios involving characteristic classes of qubit-pair 
states were discussed; in one of these, the situation in which collective dephasing, that is, dephasing between two- 
qubit-basis states occurs under a non-zero global magnetic field alone, disentanglement was not explicitly treated. 
Here, we discuss disentanglement in this case. The general conclusion in (Iol | was that bipartite disentanglement can 
be much more rapid than local decoherence under local dephasing. Here, we demonstrate a more general and precise 
result that accords with that conclusion: the decoherence timescale bounds the disentanglement timescale from above 
for both local and collective phase-damping noise environments for the two classes of states introduced there. 



A. MODEL AND MEASURES 



There exist many related two-component few-state models corresponding to physical systems the coherence of 
which may be influenced. For example, [l(| modeled two-qubit systems such as pairs of atoms in a dephasing-noise 
environment, whereas (l3j modeled two-qutrit systems corresponding to two three-level atoms whose two uncoupled 
excited levels undergo spontaneous emission to the ground state. In this section, we examine the model first introduced 
in [To| . namely, that of a two-qubit system subject to classical dephasing noise. 

Consider a two-qubit system coupled to external sources of noise that can act both on single qubits individually 
and on the joint two-qubit system as a whole. There are many physical situations that may be described in this 
way, for example, two spin-1/2 fermions on a solid-state matrix subject to random external electromagnetic fields. 
In this section, we specifically consider situations wherein the interaction between this compound system and the 
environment is described by the following Hamiltonian. 

H (t) = [B (t) (a z A + of) + b A (i) at + 6 B (t) of] , (1) 

where, for example, fi is taken to be the gyromagnetic ratio of the particle and B (t), 6a (t), and 6b (t) represent 
stochastic environmental noise from electromagnetic fields, in particular, with B(t) 0; we take fi = 1 and a z to be 
the well-known Pauli matrix; superscripts A(B) indicate the particle Hubert spaces Ha(b) hi which the operators act. 
The noise terms B (t) , 6a (i) , and 6b (t) are taken to be statistically independent classical Markov processes satisfying 

<B(i)> =0, (2) 

(B(t)B(t')) = ^S(t-t'), (3) 

(bx (*)) - 0, (4) 

(bx(t)b x (t')) = ^-6{t-t'), X = A,B, (5) 
M 

where (• • • ) is the ensemble time-average; Tx denotes the phase-damping rates associated with bx(t) (X = A, B) and 
Tab is that associated with B(t). The two-qubit standard-basis eigenstates here are 

|1)ab = | + +)ab, |2) ab = I + ->ab, 

|3)ab = I - +>ab, |4)ab = I - ->ab- (6) 

The Hamiltonian of Eq. [TJ the same one assumed in [1(3], is representative of the class of interactions giving rise 
to pure dephasing decoherence. The corresponding dynamical process is described here as a "quantum channel," in 
accord with the standard terminology of qubit evolution as described by linear maps [1, 0| , namely, a phase-damping 
channel, using the operator-sum representation which we now introduce. 

The time-dependent density matrix for the two-qubit system can be obtained by taking ensemble averages over the 
three noise fields, B (t), bA (t), and 6b (t), that is, 



p(t) = ((( Pst (*)»>, 



(7) 
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where the statistical density operator p st (t) and the unitary operator U (t) associated with H (t) are 

p at (t) = U(t) P (p)rf(t) and 



U (t) = exp 



i / dt'H{t') 



(8) 
(9) 



respectively. It is helpful to consider the dynamical evolution of p(t) as a completely positive trace preserving (CPTP) 
linear map £(p), that is, a combination of local and collective quantum channels, any of which can be turned off in 
particular cases, taking an input state p (0) to the output state p (t) given by the operator sum 



N 

p(t)=8(p(0)) = J2^lit)p(0)EAt), 
where E^ are decomposition operators that satisfy the completeness relation 



(10) 



(11) 



In each of the various cases considered here, the internal structure of the E^ will be accord with the Hamiltonian; 
various terms may or may not nontrivially contribute in a given case. The most general solution of Eq. 1 1 01 assuming 
that the system is not initially correlated with any of the three environments, is 



P® = E E (DtEfE^p(0) (E?EfD k ), 

i,j = l k=l 



(12) 



where the terms describing the interaction with the local magnetic fields &a (t) and &b (t) involve the decomposition 
operators 



Ef 



1 

7A(t) 



1 

7B (t) 



El 





w A (*) 





wb(*) 



(13) 
(14) 



By contrast, the operators describing the collective interaction with a non-zero global magnetic field, B (t), are 



Do 



7(f) 

10 

10 

7 (0 





\ u 2 (t) 







w 3 (t) 



D a = 



[Io| . The time-dependent parameters appearing in Eqs. 13-17 above are 



,7B(t) 



7aW 

7 W = e" 



-*/2T B 



t/2T AB 



wi (t) = v/1-7 2 , 

W2 (t)=- 7 yi-7 2 , 

^3 (t) = V / (l-7 2 )(l-7 4 ), 



(15) 



(16) 



(17) 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
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where Tx = ^— (X = A, B,AB) are the phase-relaxation times introduced in [lfj, Tx being the rate parameters 
appearing on tire right-hand sides of Eqs. 3 and 5; from here on, for tractability of notation, time does not explicitly 
appear as an argument for these quantities but is implied. 

The (two-qubit) entanglement is measured here by the concurrence 



C (pab) = max 0, a/Ai - \/\ 2 - \/A 3 - V A 4 , (24) 
where Xi are the eigenvalues of the matrix 

PABpAB = PAB {Oy <8> O*) p* AB g> (25) 

indexed according to decreasing magnitude, p* A „ d enotes the complex conjugate of pab, and is the standard 

Pauli matrix acting in the space of qubit A(or B) [151 ]. Concurrence is related to the canonical measure of entanglement, 
the entanglement of formation, "Ef", by 



E , w=h (i±jiE^my (26) 

where h (x) = —xlog 2 x — (1 — x) log 2 (1 — x), of which it is, therefore, a monotonic function. 

B. DECOHERENCE AND DISENTANGLEMENT 

In the standard-basis representation of Eq. 6, the generic pure state of the two-qubit system is 

I*) = o|l> AB + H2)ab + c|3) AB + d|4) A B, (27) 

a normalized state-vector with a,b,c,d € C. The generic class of two-qubit pure states represented by contains 
two pertinent subclasses, distinguished by the large-time behavior of their coherence under collective dephasing noise, 
that is, dephasing in which each qubit interacts with the same global noise field B(t): one class is fragile, whereas the 
other is robust [101 ] . 

(i) The fragile class has the forms 

\<h) =a|l)AB + &|2) A B + d|4) A B , (28) 
\<h) =a|l) A B+c|3) A B+d|4)AB • (29) 



(ii) The robust class has the forms 



=a|l>AB + 6|2>AB + c|3) A B , (30) 
\ip 2 ) =b\2) AB + c\3) AB + d\4) AB . (31) 



Without loss of generality, we consider the first of the two forms in each case, for the purposes of precisely characterizing 
the decoherence and disentanglement behavior of these classes. 

The fragile states have been shown to decohere under collective dephasing noise in such a way that all off-diagonal 
matrix elements go to zero exponentially as a function of time. A well known example of the fragile class is that with 
a = d = 1/V2, b = c = 0, namely, the spin-triplet Bell state |$ + ). By contrast, the robust states remain partially 
coherent under collective noise, some off-diagonal elements being constants [17j . A well known example of the robust 
class is that with a = d ~ 0, b — — c — l/\/2, the spin-singlet Bell state |\&~), which is known to be decoherence-free 
under collective dephasing [la]. In fact, each class of states contains two Bell states. 

In the evolution of state entanglement under local dephasing noise for both the fragile and robust classes 
was characterized; it was shown for both classes, for both the case of single-qubit local and the case of multi-local 
dephasing noise, that the disentanglement rate can be much faster than the decoherence rate when either (and so 
both) occur. We now generalize and sharpen this result, by pointing out that the decoherence timescale bounds the 
disentanglement timescale from above and explicitly showing that such a bounding relation also holds under collective 
dephasing noise, i.e., when B (t) ^ 0, a case in which disentanglement was not explicitly investigated in [ic| and only 
to a limited extent and in a somewhat different manner and context in the preceding study [9j . Here, we consider all 
possible cases in an identical way in a common formalism (that in subsequent sections is extended to the three-qubit 
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context). In this way, we confirm and sharpen the general conclusions for two-qubit systems drawn in |10l | before 
proceeding on to consider the situation of central interest here, namely, three-qubits subject to dephasing noise. 
The initial-time density matrix for the representative fragile-state class is 



Pab (t = 0) = P(\4>t)) = 



( M 

b a 


\ da 



2 a b* ad* 
\b\ 2 bd* 



whereas for the robust class one has 



Pab (t = 0) = PQih)) = 



(\a? 
ba* 
c a 






db* 



ab* 

\b\ 2 
cb* 






o \d\ 2 



a c 
be* 



(32) 



(33) 







P(\v}) indicating the projector corresponding to the generic state- vector argument \v) [l0 |. 

The time-evolved density matrix of a system beginning at t = in the fragile-class state under the collective 
dephasing channel described by Eqs. fTSHTTl is given by p{t) — E (p (0)) = J2k=i ®\p (9) ^ k '- 



I \a\ 2 a 6*7 ad*j 4 \ 



Pab (t) 



a 

b a*7 


\ da* 7 4 



\b\ 2 


d6*7 



bd*-y 


2 



(34) 



IC J 



which exhibits two pertinent two-qubit decoherence timescales, parameterized by 7 4 (fast) and 7 (slow), and becomes 
diagonal in the limit of large times. By comparison with the result in the independent multi- local noise case (see Eq. 
64 of 0) wherein two timescales of decoherence exist that are determined by 7a(b) an d 7a7b, the decay of some 
off-diagonal terms in this case is more rapid, being determined by 7 4 . 
The robust states p^ B (0) under collective dephasing noise become 



Pab it) 




(35) 



and, by contrast, do not become diagonal in the large-time limit. Furthermore, also in contrast to the case of the 
fragile class, we see that the off-diagonal elements that do decay do so according to only a single, slow timescale, 
namely, that given by a single factor of 7. The robustness of this class is due to the fact that the off-diagonal elements 
P23 and P32 (corresponding to the symmetric spin-singlet component [16]) do not decay; the robust states remain 
partially coherent in the limit of large time. By contrast, under multi-local dephasing (see Eq. 67 of [H|), density 
matrix elements P23 and p%2 do decay to zero, in the timescale determined by 7. Under collective dephasing noise, 
the two-qubit system coherence-decay timescales determined by the 7 factors in the above density matrices, are thus 

(36) 
(37) 
(38) 



F(slow) 
T 2-dcc 


-»( 


' 1 ' 

\ Tab , 


F(fast) 






T 2-dcc 


-\ 


VTab 


T 2-dec 


-i 


' 1 ~ 

,r A B, 



(for decaying elements). 

Now consider single-qubit decoherence under collective dephasing noise, corresponding to the reduction of magnitude 
of off-diagonal elements of the single-qubit reduced density matrices. Because the action of B (t) on each qubit is 
identical, it suffices to consider the reduced density matrix of just one qubit (here, we choose A) for both the fragile 
and robust classes. One finds 



Pa = tr B PAB (*) = 
Pa = tr A p A B (*) = 



a| 2 + |&r bd*<y 
db*j \d\ 2 

2 1 IJJ2 



\b\" ac'7 
ca*7 |c| 2 



(39) 
(40) 
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For the two classes, single-qubit coherence is entirely lost as parameterized by 7, namely, in the slow decoherence 
time-scale given above, so the single-qubit decoherence timescales are the same: 

(41) 
(42) 

Quantum state disentanglement is described by a reduction of concurrence, C(p). For fragile states, one finds 

C (PAB (*)) = 2 T 4 kl |d| . (43) 
Disentanglement of these states is therefore characterized by the time 

that is, the fast time-scale for two-qubit decoherence above, cf. Eq. 37. A factor of 7 4 appears here, by contrast to 
the case of multi- local dephasing (see Eq. 66 of [lj|) where disentanglement takes place via the factor 7a 7b- Under 
collective dephasing, disentanglement thus occurs more quickly than in multi-local dephasing due to an additional 
factor of 7 2 in off-diagonal terms (assuming jab, J A and jb to be comparable in value, which enables a meaningful 
comparison of effects). For the robust class of states, the entanglement is 

C(pl B (t))=2\b\\c\ ; (45) 

that is, there is no disentanglement for the robust class, which also does not decohere. We therefore see that this class 
of state can be identified on the basis of the robustness of its entanglement, as well as that of its coherence. In the 
following section, we find a similar behavior distinguishing the two major entanglement classes of three-qubit states. 

Thus, the relation between the timescales of decoherence and disentanglement described in continues to hold 
under the collective dephasing channel, in which B (t) ^ and bx(t) = 0. With the above results, together with 
those of [10|, one now can consider the behavior of entanglement under dephasing noise to be fully characterized for 
both these complementary classes of two-qubit system states, because both local-noise and collective-noise cases have 
been fully analyzed. In the following section, we consider three-qubit composite systems under local and collective 
dephasing noise, which involves a greater number of scenarios involving combinations of these noise types. 





III. THREE-QUBIT SYSTEM 

Now consider the effect of similar dephasing noise on states of three qubits. In three-qubit systems (ABC) there 
are known to be two qualitatively different pure-state entanglement classes, 

|W S ) = ai|001) + a 2 |010) +a 4 |100) , 
|GHZ S ) = a |000) + a 7 |lll), (46) 

2 2 2 2 2 

with coefficients € C such that |ai| + 1 0-2 1 + |a4 1 = 1, and | ao I + 1 0,7 1 = 1, defined by interconvertibility under 
stochastic local operations and classical communication (SLOCC) operations [la ]. We will accordingly consider each 
class individually under the available combinations of local and collective dephasing-noise that this system might 
encounter, assuming that a single noise type affects each qubit. In the above classification of (here, initial) states, the 
distinguishing property of the W class is, equivalently, the robustness of its internal entanglement under the loss of 
any one qubit of the system, in the sense of retaining the entanglement of the remaining two qubits. By contrast, the 
GHZ class has maximum three-qubit entanglement [19| but no two-qubit entanglement after the loss of any qubit. 



A. MODEL AND MEASURES 



We now consider the three-qubit system ABC in dephasing-noise environments that may act at the level of each 
of the three individual single-qubit subsystems (X = A, B, C) and/or the three two-qubit subsystems (XY = AB, 
AC, and BC), or the full three-qubit system (ABC) alone. Thus, the noise terms considered here are those due to the 
local-dephasing environment which acts on a single qubit, B", a two-qubit collective-dephasing environment, B^ 2 \ 
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and a three-qubit collective-dephasing environment, B^ 3 \ We consider only cases in which a given qubit is influenced 
by a single type of noise, in order to clearly distinguish the influence of combinations of dephasing noise on three-qubit 
systems. 

The corresponding three-qubit interaction Hamiltonian is 



2 

(2) 



+ (t) (af + 
+ Bfl (t) (at + a z c ) 
+ (t) K + ^ + .C) 



(47) 



where h is taken to be unity, and for example, \i is the gyromagnctic ratio and the B^> (i = 1,2,3) are the imposed 
stochastic magnetic fields at each of the three available scales, with subscripts indicating the region of influence in 
cases in which only subsystems are affected; Pauli matrices are also labeled by the subsystem on which they act. The 
standard three-qubit Hilbert space basis assumed here is 

|1)abc = |000)abc, |2)abc = |001)abc , 
|3)abc = |010)abc, |4)abc = |011)abc , 
|5)abc - |100)abc, |6)abc = |101)abc , 

|7)abc = |H0)abc, |8)abc = |1H>abc , (48) 

where \ijk) abc = I*) a ® |j)b ® \k)c (h j, fe = 0, 1) denote the eigenstates of the product operator <g> af £g> . 
For simplicity, much as before, we assume that the stochastic fields are taken to be classical and characterized as 
statistically independent Markov processes satisfying 



(B { x ] (*)) = 0, (49) 
(B^(t)B^(t')) = ^5(t-t'), (50) 
(B<&(t))=0, (51) 
(Bf Y (t)Bf Y {t')) = ^6 {t-t'), (52) 
(S (3) (*)) = 0, (53) 

(54) 



where (• • • ) stands for ensemble average. Here Tj is the damping rate due, for example, to the respective interaction 
with magnetic fields localized to i — 1,2,3 qubits at a time. The imposed white-noise conditions on the three-qubit 
system and any subsystems ensure a Markovian time-evolution. They also require interactions at a given scale to be 
of the same strength, which allows one to make an objective assessment of the importance of interactions by scale, as 
is done below. 

The compound-system density matrix can be obtained by taking the ensemble averages over the three noise sources, 
which are again given by decomposition-operators. For the single-qubit dephasing channel, in this section labeled T>, 
we have the following decomposition-operator matrices for qubit A at the one-qubit level. 



1 

7A(t) 

w A (t) 



1 












!)« 




l 


1 


!)• 












'(I 


1 



(55) 
(56) 
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(For the dephasing of qubits B and C, the matrices are analogous. Note also that the choice of decomposition operator 
labels in this section differs from that of the previous section, which there was made to agree with the notation of 
[Tot-) A s before, the time parameter in 7s and ujs will be implicit from here on. 

For the collective dephasing channel that acts only on two-qubit subsystems, in this section labeled £ , we have the 
following decomposition-operator matrices for the AB qubit system at the two-qubit level, which are just those that 
are collective at the two-qubit level. (For collcctive-dephasing-noise contributions acting the two-qubit systems BC 
and AC, the decomposition operator matrices are analogous.) 



/Tab 

10 

10 

V 7 AB 

/CJAB1 





V ^ AB2 



£ 3 AB 



/o 



V C^ABS 









1 
1 



1 
1 



1 
1 



(57) 



(58) 



(59) 



Finally, for the three-qubit collective dephasing channel, we have the following decomposition-operator matrices acting 
on the whole three-qubit system. 



tiABC 



/7abc 
1 
10 
1 



V 



tiABC 



( WABC1 













1 

1 























(60) 



pABC 
^3 



/o 













\ 



7ABC / 









WABC2 / 











(61) 



(62) 



WABC3 / 



The parameters appearing in the previous matrices are: 



-t/2Ti 



= y/l-e-t/T* = \/l 



Ui2 



<*3 = ^ (I -if) (I -if)- 



(63) 
(64) 

(65) 

(66) 

(67) 
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Here, as before, Tj = ji- is the phase-relaxation time associated with the corresponding interaction. Below, the 
numeric labels i may be replaced by the labels of the pertinent qubits A, AB, ABC, and so on, for specificity. 

The explicit form of the time-evolved density matrix subject to the available one-qubit, two-qubit, and three-qubit 
dephasing noise sources is given by 

P(t)= E £ (f^E^E^E^D^D^D^) p (Q) (^A^B^C^AB^BC^AC^ABc) (6g) 
i,j,k—l /,m,n,p— 1 

For each sort of dephasing-noise environment considered, the relevant off-diagonal elements from the single-qubit- 
reduced, two-qubit-reduced, and three-qubit density matrices exhibit any one-qubit, two-qubit, and three-qubit deco- 
herence effects, respectively, and bipartite disentanglement occurs when there is a reduction of (3%™, defined in Eq. 
|2"H where X, Y — A, B, C. (See the note [2(| for a discussion of the grounds for this choice of entanglement measure.) 



B. BEHAVIOR OF W-CLASS STATES 



Here, the environmental noise channel inducing single-qubit local-dephasing on qubit A is given explicitly; the 
single-qubit local dephasing noise acting on qubit B and qubit C, respectively, are of the same general form. In the 
case of subsystem-collective dephasing noise, considered in Subsection 2 below, where the collective-noise channel 
acting on qubits A and B is discussed, a similar statement applies for the channels acting on B and C only and A and 
C only. The initial generalized W-class state density matrix 



Pabc (0) = P(|W*» = 











<i2a* 



a4a| 





\ 



a\a 2 

o 

0,40,2 









a\~a\ 
a 2 al 




1 04 1 







\ 













J 



(69) 



We now examine the behavior of the time-evolved state pit), given by Eq. [|55] above, in detail. 



1. ONE-QUBIT LOCAL DEPHASING CHANNEL: V 



Under local dephasing noise, 



Pabc (*) = 











0,20,1 





^ 

o^Ija 



a\a 2 

\a 2 \ 2 a 2 al 7A 


















a 4 a^7A «4ffl27A |a.4 






I "4 















/ 



is the time-evolved three-qubit density matrix obtained when beginning at t = with a pure W-class state, 
sees that several but not all off-diagonal elements decay under local dephasing noise, as determined by 7a- 
decoherence timescale for these terms is thus 



w 

r 3-dec,I> 



= 2 



(70) 



One 
The 



(71) 



where the notation 



c(r) 
'i— dec,C> 

c(r) 
T i-dis,C 



(72) 
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represent the decoherence and disentanglement timescales, respectively, where i is the number of qubits affected, the 
pertinent dephasing noise channel is C and the state class is c; further special cases may be designated by relative 
decay rate, r ( = fast or slow), and is used from here on. The two-qubit reduced density matrices are then 



Pic (*) = 



PBc(t) 



|cii| 2 

|a 2 | 2 a 2 a|7A 

aia^jA |a4 1 2 



|a 2 | 2 

|ai| 2 aia|7A 

a 4 a*7A |a 4 | 2 

0. 

|a 4 | 2 

|ai| 2 aiCL* 2 

a 2 a* |a 2 | 2 





(73) 



(74) 



(75) 



Thus, when part of a triad beginning in a W-class state, qubit pairs that include the influenced qubit decohere 
fully, as exhibited by the presence of the factors 7a in the expected off-diagonal terms of the first two of the above 
density matrices, whereas subsystem BC maintains its coherence; no such factors appear in their joint reduced density 
matrix because in this case the dephasing channel only acts on qubit A. The two-qubit system decoherence timescale 
determined by the factors in the above density matrix is thus 



T 2 -dec,2? — L 



1 



Finally, the singlc-qubit reduced density matrices are 

Pt (*) = I 
f% (t) = 
P$ (t) = 



2 _i_ U_ |2 q 



a 2 





+ 




|a 4 

ai| 2 + M 2 



2 > 







because the individual-qubit matrices are already diagonal at time t = and so cannot be dephased further. 
Consider now the degrees of entanglement of these two-qubit subsystems as quantified by the concurrence. 



AB 



AC 



4|a 2 | 2 |a 4 | 2 7 A , 
4|ai| 2 |a 4 | 2 7 A , 



C 2 = 4|a 1 | 2 |a 2 | 2 



(76) 

(77) 
(78) 
(79) 



(80) 
(81) 
(82) 



The disentanglement timescale determined by the gamma factor for the subsystems AB and AC that decohere is 
thus 



w 

r dis,X> 



(83) 



One sees that disentanglement in subsystems containing a qubit on which local dephasing noise acts is faster than 
both the two-qubit and three-qubit decoherence: 



T dis,r> T 3-dcc,X> ) 

T dis,-D ^ r 2 -dec,X> ■ 



(84) 



Neither decoherence nor disentanglement occurs in subsystems not containing qubit A. 
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2. TWO-QUBIT COLLECTIVE DEPHASING CHANNEL: £ 

In the presence of two-qubit collective dephasing noise, acting on the two-qubit subsystem AB, the time-evolved 
full three-qubit density matrix is 



pXbc (*) 


















{) 










|3l| 2 


5i3j7AB 





01047AB 














O20J7AB 


M 2 











































4 0i7AB 


a4a 2 





M 2 



































{) 












































0/ 



(85) 



There is a decay of some but not all off-diagonal density matrix elements under the local dephasing noise, as determined 
by 7ab- A different set of terms now decay compared to the case of local dephasing noise exhibited in Eq. 70, although 
in both cases the terms P25 and decay. The timescale for the decoherence that occurs in this three-qubit system, 
determined by the two-qubit gamma factor 7ab in the above density matrix, is therefore 



r 3-dcc,£ — / 



(86) 



The two-qubit reduced density matrices obtained from Eq. 70 are 



PSc(t) = 






v 

/ M 2 



V 









V 





0,40,2 \a4\ 2 










lail 2 
a4fli7AB 




_ l_«i| 2 

02ai7AB 







5ia27AB 

|a 4 | 2 




aia^AB 

|a 2 | 2 





(87) 



(88) 



(89) 



It is significant that the off-diagonal elements of p§P do not decay, whereas the those of pffi and pyp do. By 
comparing these expressions with those of Eqs. 73-75, we also see that collective noise acting on subsystems including 
only one affected qubit has an effect analogous to that of the local noise environment on that channel. This behavior 
also appears below in the case of total-system collective dephasing of initially W-class states, where a group of qubits 
subjected to the fully collective noise also do not lose decoherence as a result of that collective noise. The timescale 
of two-qubit system decoherence determined by the factors in the above density matrices, when it occurs, is thus 



r 2-dcc,£ I 



(90) 



Again single-qubit reduced density matrices are diagonal from the outset, cf. Eqs. 77-79, and so cannot decohere 
further. 

The degrees of entanglement of two-qubit subsystems are 



C 2 



AB 



r 2 



y AC 

Cl c = 4|oi 



4 |a 2 | 2 |azi| 2 , 
4|ai| 2 |a 4 | 2 7AB 



|02| 7ab • 



(91) 
(92) 
(93) 



Comparing these expressions with the local-dcphasing noise scenario expressions in Eqs. 80-82, one again sees sim- 
ilar behavior when only one qubit of a subsystem is subject to collective noise. The timescale of disentanglement 
determined by the factors in the above density matrix, when it occurs, is thus 



(94) 
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In this case, one again sees that disentanglement, which occurs only when decoherence occurs, is always faster than 
three-qubit and two-qubit decoherence in all instances in which these effects occur: 



T dis,£ ^ r 3-dec,£ 



W 
r dis,£ 



< T, 



W 

2-dec,£ 



(95) 



3. THREE-QUBIT COLLECTIVE DEPHASING CHANNEL: T 

The time-evolved three-qubit density matrix, however, shows no effect from the three-qubit collective dephasing 
channel: 



Pabc (*) = 



/ 

o M 2 

a,2a\ 


04 a \ 



V 



0\ 

aiaj 
a 2 a% 



ai<3,2 
|a 2 | 2 









|a 4 | 














000/ 



(96) 



Hence, the W-class states can be identified as a class that is entirely robust under fully collective noise, in the sense 
being "decoherence-free," which can be attributed to the permutation symmetry introduced with the initial state 
IVF 9 ), analogously to what occurs in two-qubit systems under collective dephasing (cf. [IB]). For the GHZ-class 
states, however, one finds decoherence does occur, as shown in Section C, below. 



4. THREE-QUBIT MULTI-LOCAL DEPHASING: V A V B V C 



Under multi-local dephasing noise, wherein each qubit is individually subject to the same sort of noise, one finds 
the time-evolved density matrix arising from initially W-class states to be 



pXbc (*) = 



/ 











o\ 


|ai| 2 


aia27B7c 


aia|7A7c 


000 


a 2 a*7B7c 


l« 2 | 2 





a 2 a|7A7B 


000 














000 


a 4 a^7A7c 


04a27A7B 





M 2 


000 














000 














000 


V 











000/ 



(97) 



wherein all off-diagonal elements decay according to a product of two differing 7 factors. The three-qubit system 
decoherence timescale determined by those factors (recalling that 7^ = 73 = 7c) in the above density matrices is 
thus 



w 



T 3-dec,X' A 'D B I5 c 



(98) 
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The two-qubit reduced density matrices are 



Pbc (*) 



/\ai\ 2 













|a 2 | 2 


a2047A7B 


° 




a 4 a27A7B 


|a 4 | 2 




\ (J 


U 







u / 


/ |a 2 | 2 








°\ 





lail 2 


aia47A7c 








a4«i7A7c 


|a 4 | 2 





V o 








0/ 


/ l«4| 2 








°^ 





l«i| 2 


aia27B7c 








«2at7B7c 


|a 2 | 2 





V o 








0/ 



(99) 



(100) 



(101) 



similar in form to those in Eqs. 73-75 but having a 7x factor in off-diagonal elements for each corresponding local (to 
qubit X) noise source. One sees that all off-diagonal elements decay, as determined by two pertinent 7 factors. The 
two-qubit system decoherence timescale determined by the factors in the above density matrices is similarly thus 



w 

T 2-dec,V A V B V c 



(102) 



The single-qubit reduced density matrices are again diagonal, as shown in Eqs. 77-79, and so are not able to dephasc 
further. 

The degrees of entanglement of the two-qubit subsystems given by concurrence are 



C 2 

c 2 



us = 4|a 2 | 2 M 2 7a7b > 
u = 4|ai| 2 |a 4 | 2 7 A 7c » 



ril a I- 1 2 1- 1 2 2 2 

Cbc = 4 l«l I M 7b7c . 

with the particular 7 changes reflecting the dephasing channel. The disentanglement timescale is thus 



(103) 
(104) 
(105) 



T 2-&is,V A V B V c " 



1 



2 Vri 

We thus find again that disentanglement is always faster than three-qubit decoherence and two-qubit decoherence 



(106) 



T dis,I5 A I5 B r> c ^ T 3-dcc,-D A X> B £ c 

_W . W 

T dis,-D A 2? B 2? C ^ T 2-doc,r> A D B X> c 



(107) 



for multi-local dephasing noise. 



5. ONE-QUBIT LOCAL, TWO-QUBLT COLLECTIVE DEPHASING T>£ 

For three-qubit systems, the possibility of combinations of local and subsystem-collective dephasing noise exists, 
which we now consider. Under such noise, where qubit A, for example, is affected by local noise and the remaining 
subsystem BC is subject to collective noise (i.e. noise collective at the scale of BC), the time-evolved state is 



PABC (*) 

















\ 





l«ll 2 


aic 2 





aiC47A7BC 








a 2 ci 


|a 2 | 2 





fl2C47A7BC 


























a4Ci7A7BC 


a4C 2 7A7BC 





N 2 









































\o 














/ 



(108) 
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The effective three-qubit system decoherence timescale determined by the factors in the above density matrices is 
thus 

^3 W dcc^ - +2^-) . (109) 

In this case, one finds slightly different decay behavior than in cases above: the total decay rate is composed of a 
product of 7 factors. In this case and in subsequent similar cases, we will assume, for clarity of exposition, that cither 
one of the two 7 factors, one corresponding to local noise and the other to collective noise, is much larger or that the 
two are of comparable strength. 

The two-qubit reduced density matrices in this case are 

|ai| 2 

Pab (*) = I n - ~' a2? ~ a2 VT C n I' < ll( " 
a 4 c 2 7A7BC M 1 



a 2 | 2 

P f c (t)= I 9 _ J Sl ' 2 1 ni) 

! &4C17A7BC \ai\ 














a 4 | 2 














\ai\ 2 


a\c 2 








a 2 ci 


l« 2 | 2 


















*("-! 5 S£ ffih (112) 

V 00/ 

The off-diagonal elements also decay due to the appearance of a combination of the two pertinent 7 factors in the 
first two subsystems above. Because the two qubits of subsystem BC are subjected to collective dephasing noise 
and because W-class states are symmetrized, this subsystem does not decohere, differing from the case of two-qubit 
subsystem-collective dephasing only and the three-qubit collective dephasing cases. The two-qubit system decoherence 
timescale, in subsystems AB and AC, determined by these factors is thus 



r?- Acc ,v £ =\- 1 ) +2 {r~ 2 ) ' (113) 

the same combination of rates as for three-qubit decoherence. 

Recall again that the single-qubit reduced density matrices, given by Eqs. 77-79, are fully dephased from the outset. 
The entanglement of subsystems given by the concurrences 

C AB = 4|a 2 | 2 |a 4 | 2 7A7lc , (114) 

Ci c = 4|a 1 | 2 |a 4 | 2 7A7BC . (H5) 

Cl c = 4|a 1 | 2 |a 2 | 2 , (116) 

has behavior determined by the factors in the above density matrices, so that for subsystems AB and AC 



Disentanglement is therefore always faster than two-qubit and three-qubit decoherence under this noise, 

T dis,V£ ^ r 3-dcc,I5£- i 
T dis,V£ < T 2-dec,VSi K LLii ) 

when it occurs. 



6. DISCUSSION 



From the above results for the W class of states, we have shown that the disentanglement timescale is always faster 
than or equal to the decoherence rate as described by the off-diagonal terms in all five dephasing-noise environments, 
namely, one-qubit local (V), two-qubit collective {£), three-qubit collective (J 7 ), three-qubit multi- local (DVD), and 
one-qubit local plus two-qubit collective (D£ ) environments. Thus, the decoherence rate bounds the disentanglement 
rate from below for systems beginning out in W-class states. 
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C. BEHAVIOR OF GHZ-CLASS STATES 

The analysis of the GHZ-class of states will proceed in a similar fashion as that of the W-class but can be explicated 
more briefly. The initial density matrix of the GHZ-class state of Eq. 46 is 



7i.\2 



P^(0) = P(|GHZ*» 



a a 7 \ 













V a 7 a* |a | 2 / 



/ ho 









(119) 



1. DEPHASING CHANNELS 

The effect of all dephasing noise types upon the GHZ-class density matrix is similar, differing only by the represen- 
tative decay factor 7 of the decay channel, in the five previously considered noise-scenarios: 



V A : 


7 


£ AB : 


7 


T : 


7 


A T> B V C : 


7 


v a £ bc . 


7 



1 — qubit local, 
2 — qubit local, 
3 — qubit collective, 
7a7b7c multi — local, 

1 — qubit local and two • 



4 

7ab 
7abc 



7a7bc 



qubit collective 



(120) 



which appear in all of the off-diagonal elements in the density matrices in the following time-evolved density matrix, 
so that 



RtTc (*) = 



I |a | 2 a d^\ 

000000 













V a 7 a^ 








|a T | 



(121) 



/ 



which decays to a diagonal matrix as determined by the particular form of 7 under the influence of the noise models 
introduced in the previous section. The three-qubit-system decoherence timescales are thus 



GHZ 
T 3-dcc.I? A 



GHZ 

T 3-dec,£ AB 



GHZ 

T 3-dec,^ AB C 



GHZ 
r 3-dec,X> A X> B X> c 



GHZ 

T 3-dec,X> A £ B C 



1 

2 
1 

2 
2 
3 

= 2 



1 

1 
1 

rT 
1 



(122) 
(123) 
(124) 
(125) 
(126) 
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The two-qubit reduced density matrices are, unlike the case of the W class, seen to be diagonal from the outset and 
so cannot decohere or dephase further: 



P^ Z (t) = 



PIT (t) = 



ao\ 2 
















































I — I 9 

H 2 


/ 


ao| 2 











\ 














) 

























l«r| 2 




N 2 







































) 











larl 2 





The single-qubit reduced density matrices are similarly always diagonal: 

12 o 



pT z (t) - 
pT z (t) - 

pg HZ (t) 

For two-qubit-subsystem entanglement, one finds 



I CIO 





\ a ?\ 
2 



|a 7 P 
2 



\a 




\ a 7\ 



r 2 

^AB 


= o , 


r 2 

°AC 


= o , 


r 2 

°BC 


= . 



(127) 



(128) 



(129) 



(130) 
(131) 
(132) 



(133) 
(134) 
(135) 

The reduced states of all proper subsystems in states initially of the GHZ-class are neither coherent nor entangled 
at any time, when subject to all the above forms of dephasing noise; all entanglement appears only at the triple-qubit 
level, for which there is no analytic mixed-state entanglement measure of which we are aware. Nonetheless, one 
sees that, at the two-qubit level, the behavior of this class of states still formally agrees with our general conclusion 
regarding the W-class, with the difference that the two-qubit (lack of) decoherence trivially bounds the trivial (lack 
of) bipartite entanglement. Unfortunately, as in the case of the W-class states but more significantly in this case, the 
lack of a well defined analytic mixed-state three-qubit entanglement measure precludes one from performing a detailed 
comparison of three-qubit disentanglement versus decoherence at the three-qubit level, which would be useful. 



IV. CONCLUSIONS 



We have shown for the fundamental three-qubit entanglement classes that there are different timescales for two- 
qubit disentanglement and decoherence under dephasing noise. In particular, we have shown for three qubits initially 
in an entangled state subject to an array of combinations of dcphasing-noisc types, disentanglement, when it occurs, 
does so on a timescale shorter than or equal to that of decoherence. In the course of obtaining these results, we 
determined the precise dephasing and bipartite disentanglement behavior for each of the various pertinent classes 
of two-qubit and three-qubit states under all dephasing noise scenarios available at their respective scales, and have 
clearly identified dephasing-free and disentanglement-free classes. The coherence and bipartite entanglement behavior 
for these classes has thereby been provided for future investigators, who can now select and/or engineer states for 
practical applications, such as quantum metrology, communication and information processing, according to the 
coherence and entanglement characteristics these tasks require. 
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